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A rather general enhancement of superconductivity is demonstrated in a hybrid structure consist- 
ing of submicron superconducting (SC) sample combined with an in-plane ferromagnet (FM). The 
superconducting state resists much higher applied magnetic fields for both perpendicular polarities, 
as applied field is screened by the FM. In addition, FM induces (in the perpendicular direction to 
its moment) two opposite current-flows in the SC plane, under and aside the magnet, respectively. 
Due to the compensation effects, superconductivity persists up to higher applied currents. With 
increasing current, the sample undergoes SC- "resistive" -normal state transitions through a mixture 
of vortex-antivortex and phase-slip phenomena. 

PACS numbers: 74.78.Na, 74.25.Dw, 74.25.Op, 74.25. Sv. 



Over the last years, the superconductor-ferromagnet 
(SC-FM) hybrids received a lot of attention as one of 
the rare systems where ferromagnetism and singlet su- 
perconductivity coexist (for review, see 0). These hy- 
brid structures are looked upon as candidates for futuris- 
tic nanoelectronics, combining superconducting circuits 
with magnetic storage elements. As better understand- 
ing is needed, the ongoing studies are mainly focused on 
fundamental properties of nanoscale SC-FM samples and 
plethora of related phenomena. 

For example, although ferromagnetism in general sup- 
presses superconductivity, direct SC-FM coupling ap- 
pears to be crucial for the 7r-phase state with the critical 
current inversion in SC-FM-SC junctions and Joseph- 
son current enhancement in SC-FM tunnel structures 
with very thin FM layers 3] . On the other hand, the non- 
trivial interplay between magnetism and superconductiv- 
ity can be achieved even if SC and FM are not electron- 
ically coupled, as they still interact through the emerg- 
ing magnetic fields. In that respect, arrays of submicron 
magnetic particles are used for applying well-defined lo- 
cal magnetic fields in the underlying superconductor 
One of the first applications of these nano-magnets was 
to engineer the pinning force of superconducting films, 
such that the critical current function of applied 

magnetic field is increased due to a collective locking of 
the flux lattice to the magnetic array 0, 0|- Since then, 
because of the technological relevance, enhancement of 
critical parameters in SC-FM heterostructures is of vast 
theoretical and experimental interest. Genenko et al. 
predicted theoretically an increased edge barrier critical 
current in superconductors completely surfaced by mag- 
netic material . In that demagnetized magnetic 
layer acts as a magnetic screen, effectively shielding the 
Meissner state. Two years ago, Lange et al. measured 
higher critical field in SC films regularly structured by 
out-of-plane magnetized dots Q. However, this behav- 
ior strongly depended on the polarity of the applied field 
Hext'- for given FM- magnetization M, an enhancement 



of the critical parallel field (H ex t || M) was achieved 
at expense of the antiparallel one. The same behavior 
was found both experimentally || and theoretically 
in mesoscopic SC disks with out-of-plane FM dot on top. 

The first objective of the present Letter is to design 
a SC-FM hybrid structure where most critical proper- 
ties can be tailored practically at will. For that matter, 
we consider a thin submicron superconducting sample 
with a ferromagnetic dot with in-plane magnetization on 
top (see Fig. QJ. Such a device realization offers full 
exploitation of the magnetic flux pinning [lfijl. dy nam- 
ical properties of mesoscopic superconductors [llj, and 
related vortex-antivortex phenomena Due to the 

opposite magnetic field at the poles of the magnet, the 
field-compensation effects lead to the critical field en- 
hancement for both positive and negative applied field 
(see Fig. H). At the same time, our SC-FM sample might 
act as a current compensator as well: as a novel concept, 
the applied current is met by opposing FM-induced su- 
percurrents, resulting in a larger critical current. 

In our theoretical treatment of this system, we rely 
upon the Ginzburg-Landau (GL) formalism. In, the sta- 




FIG. I: A superconducting (SC) sample underneath a square 
ferromagnetic (FM) dot with in-plane magnetization M (sep- 
arated by an insulating layer). Depicted directions of M, the 
applied external magnetic field H ex t and current j ea , t are de- 
noted as positive throughout the article. 



tionary case, we solve self-consistently two GL equations, 
derived from the Gibbs energy functional. For all details 
of this approach, we refer to Refs. (l^Hsf. 

To understand the dynamical properties of the device, 
we studied the current-voltage characteristics using the 
time-dependent Ginzburg-Landau (TDGL) equation Q 

u fd. r 2 a|*| 2 \ 
= (v - iA) 2 * + (i - t - |*| 2 )*, (1) 

coupled with the equation for the electrostatic poten- 
tial Aip = div (Im(**(V - iA)*)). Here, the distance 
is measured in units of the coherence length £(0), \& is 
scaled by its value in the absence of magnetic field ipo, 
time by tgl(0) = Trh/8kBT c u, vector potential A by 
cfi/2e£(0), and the electrostatic potential by h/2eTGL(0). 
r = 2rEipo/h, with te being the inelastic electron- 
collision time. For Al samples te ~ 10ns, which results 
in r sa 1000. Parameter u — 5.79 is taken from Ref. 
[lij . Note that in Eq. 0) the screening of the magnetic 
field is neglected, as we restrict ourselves to thin SC sam- 
ples (d < £). The points where external current jext is 
injected in the sample (see Fig. were simulated as 
normal metal-superconductor contacts, i.e. with W = 
and — V</? = jext- At the remainder of the sample edges, 
Neumann boundary condition was used (j± =0). 

We consider a square Al sample with parameters eas- 
ily achievable with modern lithographic techniques: size 
o-x = o, y = 1.5/im, thickness d = 80nm, separated by an 
oxide layer of thickness I = 20nm from the square FM 
with size w x — w y = 800nm and thickness D = 50nm. 
The SC material is characterized by its coherence length 
at zero temperature, which we take as £(0) = lOOnm 
(typical value for mesoscopic Al samples j8|), and FM 
material by its saturation magnetization M. In our cal- 
culations, the FM is positioned at the center of the SC 
square, and since its stray field has opposite polarity at 
the FM-poles (see inset (a) in Fig. (2J, the total flux $fm 
penetrating the SC equals zero. This feature inevitably 
leads to the appearance of vortex-antivortex configura- 
tions for sufficiently strong magnetization M |l2^ . 

Fig. El shows the Gibbs free energy 0] of the super- 
conducting state, obtained after sweeping up/down the 
FM-magnetization, where the number of induced vortex- 
antivortex (VAV) pairs is denoted by Roman numbers. 
Note that these VAV states are the first found vortex 
states with zero total vorticity in finite mesoscopic SC 
samples. Insets (b,c) in Fig. |21 show the |^| 2 -density 
plots of successive VAV states. As vortices and antivor- 
tices are confined at the FM-poles (where the stray field 
is maximal), they are effectively kept apart by the FM. In 
other words, the superconducting region under the FM 
always remains (anti)vortex-free. As a result, supercon- 
ductivity can be sustained in the sample up to very large 
FM-magnetization (as the slope dQ / dM decreases in Fig. 
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FIG. 2: The Gibbs free energy diagram (Qo = H 2 V/8ty) as a 
function of FM-magnetization. Roman numbers denote num- 
ber of FM-induced vortex-antivortex pairs. Inset (a) illus- 
trates the FM-stray field lines; (b) the Cooper-pair density 
[upper figure, darkest color - zero density (online blue/red- 
low/high density)] and superconducting phase contourplot 
[gradation of grey color shows the circulation of phase — 2n 
(online blue/red-0/27r phase)] for state I; (c) the \ip\ 2 -density 
plots (scale adjusted for clarity) under the positive pole of the 
magnet, for states II to V (FM-edge depicted by white lines). 

|2J ■ Note that this is not the case if FM has perpendicular 
magnetization, when total flux <&fm captured by SC is 
positive and FM-induced vortices destroy superconduc- 
tivity in the heart of the sample The experiment of 
Ref. [8j revealed that when such a sample is exposed to 
homogeneous external field H ext , the H ext —T boundary 
is shifted towards positive fields due to the compensation 
with $fm, resulting in higher positive critical field (and 
consequently reduced negative one). 

In Fig. |3| the H ext — T superconducting/normal (S/N) 
phase boundary of our sample is shown, in the case of 
FM with bulk Co magnetization M = 1400G (solid line), 
compared to the case without FM (dashed line). The 
S/N boundary exhibits three novel features: (i) the M.- 
shaped boundary - the critical temperature for H ext — 
is reduced (T cm o)> an d is maximal for two symmetric 
non-zero H ext values (T cm ); (ii) for T cm0 < T < T cm 
the N-S-N-S-N multi-reentrant behavior is observed 
during H ext sweep; and (hi) substantial critical field 
enhancement is found for both H ext polarities. 

The physical reason for these phenomena lies in the 
magnetic field compensation. In this particular case, for 
Hext = 0, the S/N transition at T — T cm o occurs for 
2 vortex-antivortex pairs induced by the FM (state II 
in Fig. 0). Although their centers are pin-pointed at 
the FM-poles, these (anti)vortices are covering the whole 
sample as £(T) becomes large. When H ext > is ap- 
plied, the external flux is "absorbed" by the FM-induced 
antivortex. This effectively recovers superconductivity, 
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FIG. 3: The superconducting/normal (S/N) phase boundary 
for a square SC sample, without (dashed curve) and with 
(solid curve) an in-plane FM dot on top. The total vorticity 
of the sample is denoted by L, while indices a, b show the 
number of vortices and antivortices (induced and/or pinned 
by FM) at corresponding FM poles, respectively (L — a — b). 



and increases T c . Each kink in the Little-Parks- like S/N 
boundary with increasing H ext corresponds to a change 
in total vorticity of AL — 1, where external flux lines 
are first annihilated by FM-induced antivortices, and in 
the absence of antivortices pinned on the positive pole 
of the FM, where the stray field and H ext are aligned 
(as found in Refs. In the latter case, each ad- 

ditional vortex suppresses superconductivity and T c de- 
creases. However, the SC state remains protected at the 
opposite pole of the magnet which results in significantly 
higher critical field. In Fig. [21 an enhancement as high as 
~ 40% is achieved. This percentage can be even larger, 
if stronger magnetic materials are used (note high M cr 
in Fig. |2J). The reduced zero-field critical temperature 
T C m0 is only within few percent from T c q, but maximal 
T cm (at H ext — 1.97mT) is several percent higher than 
the corresponding T c value in the absence of the FM. 

For H ext < 0, the scenario is completely analogous, 
and S/N boundary is therefore symmetric. Note that 
this symmetry directly leads to feature (ii) , which is actu- 
ally a very rare magnetic field-mdwcerf-superconductivity 
(FIS) phenomenon. Similar unconventional behavior was 
reported earlier for materials like (EuSn)MogS8 and A- 
(BETS) 2 FeCl 4 ,and for SC films with out-of-plane FM- 
arrays on top [jj. However, our sample holds a unique 
property - FIS is achieved for both perpendicular polariza- 
tions of applied field H ext (e.g. Fig. 0for T = 0.985T c0 ). 

Obviously, the above described phenomena are directly 
related to the strongly inhomogeneous FM-stray field 
(with zero average). Yet another interesting feature can 
be found in the stray- field-induced currents j_fm- Due 
to the field landscape, these currents are actually cir- 
culating around the poles of the magnet, which ulti- 



mately results in two opposite current flows, under and 
aside the FM, in a direction perpendicular to the FM- 
polarization (see Fig. Ufa)). Obviously, with increas- 
ing FM-magnetization, the amplitudes of Sfm grow (as 
shown in Fig. Efb)), until the depairing current is reached 
and vortex-antivortex (VAV) pair nucleates at the FM- 
poles. Consequently, jfm completely reverses (Fig. ^b), 
yellow line), but changes polarity again with increasing 
M, before the appearance of the following VAV pair (see 
Fig. [21). This dual, step- like jfm profile may strongly 
affect the response of the device on the applied current 
in the ^-direction (see Fig. . In order to investigate the 
critical current and dynamical properties of the system, 
we employ TDGL formalism. The key results are shown 
in Fig. as differential resistance (obtained from calcu- 
lated I-V characteristics) as function of applied current 
j ex t- Two critical currents, denoted as j c \ and j C 2 in Fig. 
|SJ can be identified. j c % is the current at which the sam- 
ple loses its zero-resistance and transits to the so-called 
"resistive" state. j C 2 has the more conventional meaning 
of the current at which the SC state becomes unstable. 

When external current is applied to a plain SC square 
(Fig. Efa)), it is non-uniformly distributed in the sample, 
with its maxima at the side-edges (see cartoon in the in- 
set). It is at these weak points where the vortex nucleates 
when the depairing current is reached (for correspond- 
ing j ext = jci)- Due to the Lorentz force, this vortex 




FIG. 4: (a) Vectorplot of the SC-current jfm induced by an 
in-plane FM (Meissner state), superimposed on the contour- 
plot of jp M component, (b) The jf^M-profile in the central 
cross-section for different FM-magnetization (T = 0.97T c o)- 
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FIG. 5: Calculated differential resistance of the sample as function of the applied DC current, for different FM-magnetization 
(T = 0.97T c o). White areas denote Meissner state, light grey (yellow online) - "resistive", medium grey (orange) - vortex- 
antivortex, and dark grey (green) - normal state. Insets illustrate schematically the distribution of the applied (solid line) and 
FM-induced current (jf-M component, dashed line) across the sample. 



is expelled across the sample, and subsequently nucleat- 
ing again. To prevent the destruction of superconductiv- 
ity, the SC phase exhibits a jump of 2ir at this "phase- 
slip line" [15j . This phase-slip can be theoretically inter- 
preted as an infinitely-fast moving vortex. With contin- 
uous increase of j ext , the current further suppresses the 
order parameter at the contacts (see Fig. ^) finally estab- 
lishing a normal path between them (for j ext = j C 2) and 
normal- metal resistance is reached in Fig. Ufa). Note 
that the presence of S/N contacts in our simulation gives 
small, but finite resistance even for j ext < j c \. 

In the case of our SC-FM device, for small M, the 
FM induces current opposite to the external one, along 
the sample edges. Therefore, the resulting current at the 
weak points is decreased, and the critical current for vor- 
tex entry (phase-slip and non-zero resistance) increases. 
Fig. |SJb) shows exquisite j c \ enhancement of ~ 35%, for 
M = 100G (see Fig. E{b) for j FM profile, and Ref. 
However, for larger FM-magnetization, after superposi- 
tion of j ex t and jFM, although edge current is further 
suppressed, the depairing current may be reached under 
the FM where jpM is maximal. As a result, VAV pair nu- 
cleates, and SC transits to the resistive state for a low j c \ 
value (Fig. 0[c)). Following VAV creation, jpM changes 
polarity, and now compensates j ex t between the contacts 
instead of at the edges. Eventually, with further increase 
of jext, jFM is overwhelmed by the applied current, and 
VAV pair is expelled from the sample; the influence of FM 
becomes negligible as further scenario resembles the one 
of Fig. I3a): current is again maximal at the edge, phase- 
slip occurs, and superconductivity is destroyed. However, 
due to described VAV nucleation and current compensa- 
tion between the contacts, for M — 450G we obtained a 
remarkable enhancement of j C 2 of ~ 21.5% (Fig. EJc)). 

For higher M , FM may induce a VAV pair in the sam- 
ple (Fig. Efb)). In that case, even for very low applied 
current, a finite resistance was found {j c \ = 0). This 
feature can serve as a tool for experimental detection of 



VAV pairs in contrast to the Meissner state in SC-FM 
hybrids. For certain value of j ex t, vortex and antivortex 
are depinned and leave the sample, followed by an imme- 
diate phase-slip and consequent transition to the normal 
state. Due to the absence of the zero-resistance state, 
both jd and j c i are significantly decreased (Fig. Efd)). 

In conclusion, we proposed a SC-FM device where both 
critical field and critical current can be substantially en- 
hanced. Although our dynamical simulations are valid 
only in close vicinity of T c , the main idea is generally ap- 
plicable. Detailed influence of parameters and different 
dynamic regimes will be analyzed in a separate article. 
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